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Abstract 

The instantaneous Bethe-Salpeter equation, derived from the general Bethe-Salpeter 
formalism by assuming that the involved interaction kernel is instantaneous, represents 
the most promising framework for the description of hadrons as bound states of quarks 
from first quantum-field-theoretic principles, that is, quantum chromodynamics. Here, 
by extending a previous analysis confined to the case of bound-state constituents with 
vanishing masses, we demonstrate that the instantaneous Bethe-Salpeter equation for 
bound-state constituents with (definitely) nonvanishing masses may be converted into 
an eigenvalue problem for an explicitly — more precisely, algebraically — known matrix, 
at least, for a rather wide class of interactions between these bound-state constituents. 
The advantages of the explicit knowledge of this matrix representation are self-evident. 
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1 Introduction 

In principle, the appropriate tool for the description of bound states within relativistic 
quantum field theory is the Bethe-Salpeter formalism. However, attempts to solve the 
Bethe-Salpeter equation face several, well-known obstacles. In view of this, in practical 
applications usually the easier-to-handle instantaneous approximation for the involved 
interaction kernel is considered. 

In a recent paper Q, we introduced, for the somewhat simpler example of massless 
bound-state constituents, a technique for converting the instantaneous Bethe-Salpeter 
equation into an eigenvalue problem for an explicitly given matrix. Here, the analysis of 
Ref. |I| is extended to the case of nonvanishing masses of the bound-state constituents. 



2 Bethe-Salpeter Equation as a Matrix Equation 

Our goal is to demonstrate the possibility of converting the Bethe-Salpeter equation to 
a matrix equation. Thus, let us accept the same simplifying assumptions as in Ref. [![: 

• The propagators in the Bethe-Salpeter equation may be approximated by free 
propagators with some kind of effective masses of the bound-state constituents. 

• The bound-state constituents have equal masses. 

Paralleling the discussion in Ref. we consider fermion-antifermion bound states 
with spin J, parity P = (— 1) J+1 and charge-conjugation quantum number C = (— 1) J . 
The corresponding equal-time Bethe-Salpeter amplitude, or "Salpeter amplitude," x 
involves two independent components, and \&2- For two fermions of equal masses m 
and internal momentum k, in the bound state's rest frame it reads in momentum space 



X(k) 



*i(k) 



rn 



"y ■ k 



E(k) 



+ ^ 2 (k) 7 C 
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E{k) = \Jk 2 + m 2 , k = |k| , is the energy of a free particle of mass m and momentum k. 
We confine ourselves to the case J = 0, that is, to bound states with spin-parity-charge 
conjugation assignment J PC = 0~ + , denoted by 1 So in usual spectroscopic notation. 

By expanding the Salpeter amplitude into some convenient set of basis matrices in 
Dirac space and after factorizing off the (vector) spherical harmonics depending on the 
angular variables, the instantaneous Bethe-Salpeter equation for fermion-antifermion 
bound states can be reduced to a set of coupled equations for radial wave functions ||. 
For pure time-component Lorentz- vector interactions, i.e., 7°®7° kernels, it reads pL EI 



2E{k) ^ 2 (k) + 



dk' k 



rl 



V (k,k')^ 2 (k') = M^ 1 (k) , 



+ 



/ U'2 



dk'k 
W) 



'" v (k,k')^- + J- Vl (k,k') ?/ 



E{k) 



E(k') E(k) 



E(k') 



^i(fc') = M^ 2 (k) 



where, expressed in terms of a static interaction potential V(r) in configuration space, 



V L (k,k') = 8tt J &rr 2 V{r)i L (kr)j L (k'r) 
b 



L = 0,l,2, 



2 



Here, j n (z) (n = 0, ±1, ±2, . . .) are the spherical Bessel functions of the first kind |3J]. 

Expressing, for M^O, from the first of Eqs. (fj) the component in terms of the 
component ty 2 of the Salpeter amplitude x an d inserting this into the second of Eqs. ([1]) 
yields an eigenvalue equation for \l>2 with the bound-state mass squared as eigenvalue: 



M 2 V 2 (k) = AE 2 (k)^ 2 (k) + 2E(k) J 



1 U'2 



dk' k 



2 V (k,k')Mk') 



+ 2 



+ 



dk' k' 2 
~(2^f 

dk' k' 2 



'"' V Q {k 1 k') + ^V 1 {k 1 k') 



E(k) 



E(k) 



* 2 (&') 



E(k) 



V (k,k') 



m 



k 



E(k') E(k) 



V 1 (k,k r ) 



k' 



E(k') 



x 



n 1J12 



dk" k 



V (k',k")V 2 (k") 



In order to convert this eigenvalue equation to matrix form, we introduce suitable 
sets of basis functions for the Hilbert space L 2 (R + ) of (with weight function w(x) = x 2 ) 
square- integrable functions f(x) on the positive real line R + (see [0]). For a given value 
I — 0, 1, 2, ... of the angular momentum of their counterparts in three dimensions, the 
basis functions are called (j)\ (r) in configuration space and 4>f\p) in momentum space. 
They are the same as in Ref. [|TJ apart from the fact that the real variational parameter 
H > there is replaced here by the mass m of the bound-state constituents and that, 
consequently, normalizability of the basis vectors demands m > 0; their main features 
are summarized in Appendix Expanding ^ 2 (p) in terms of the radial basis functions 
4>f^ (p) , the instantaneous Bethe-Salpeter equation is solved by diagonalizing a matrix: 



Mi 

with the abbreviations 



Aij + Bij + C, 



(i) 



C 



(2) 



D 



(i) 



D 



(2) 



.4 



4 / dA: k 2 E 2 (k) 4 0) (k) (t>f (k) , 



dkk 2 E{k)(t)f\k) I dk'k^Voik^')^^') , 



r0) - 



D {1) = 
^ij — 



D 



(2) 



2 m 2 

w? 

2 

2 

m 

1 



dA: k 2 
Ejkj 



4°\k) I dk'k^V^k^k')^^') , 



(o), 



dA: A; 3 

W) 

dk k 2 
W) 

dkk 3 
{2n) 4 J E{k) 



4°\k) I dk'k'^V^k^k')^^' 



(0), 



^\k) 



dk' k 



1 1.12 



E{k' 



-V (k,k') I dk"k" 2 V Q {k',k")<j>f\k") , 



<t>t\k) J%^Vi(k,k') Jdk"k" 2 V (k',k")<pf\k") . 
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The solution of this eigenvalue problem (of course, only for a finite matrix size d) 
proceeds, step by step, exactly along the lines presented in much more detail in Ref. fl[]. 
First of all, we introduce the matrix elements of the square E 2 of the kinetic energy: 

oo 

Kijim) = J dkk 2 E 2 (k)c/)f ) (k)(f ) f\k) . 
o 

Furthermore, in order to evaluate the terms B^, c!>j , . . . , D^' analytically — which 
requires repeated applications of the Fourier-Bessel transformations (|6]) — , we have to 
expand several expressions involved in these integrals in terms of the appropriate set of 
momentum-space basis functions (pf\k), £ = 0, 1 (cf. Eqs. (28), (31), (32) of Ref . ||T|): 



E(k)< 


Pf\k) 


N 

= EM 

3=0 


m) 


<f>f\k) 


k 

E(k)' 




N 

= E C ji 
3=0 


<t.f 


(k) , 


ki 




N 

= E d Ji 

3=0 


(m) 


4>f\k) 


1 

E(k) ( 


i>f{k) 


N 

= E e a l 

3=0 




<Pf\k) 



As consequence of the orthonormality of the momentum-space basis functions <p\ (p), 
the expansion coefficients bij(rn), Cij, dij(m), and &ij{vri) may be expressed in the form 
(cf. Eqs. (29), (33), (34) of Ref. §) 

oo 

6y(m) = J dkk 2 E(k) (t)f\k) <j)f\k) , 






dijim) = J dkk 3 <pf\k)^f\k) , 
o 

These expansion coefficients are, of course, not independent but satisfy several (only in 
the limit N — > oo, exact) relations of the kind 

N N 

b ri (m) b rj (m) = E d* ri {m) d rj (m) + m 2 = K^m) , 

r=0 r=0 
N 

E b ri (m) e rj (m) = 8^ , 

r=0 
N 

E d ir {m) e rj (m) = . 

r=0 
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Employing these relations in order to investigate systematically the errors induced by 
the truncations of the expansion series in Aiij, one finds that, for instance, for d — 15 
(i.e., 15 x 15 matrices) and for N = 49 (i.e., a truncation to the first 50 basis vectors), 
all the above relations are satisfied with relative errors less than 3%. 

Finally, we'll need the expansions of the expressions V(r) 4>f\r) in terms of tfif (r): 

N 

V (r) # (r) = £ V$> (m) <f>f (r) , £ = 0,1. 

3=0 

Quite obviously, here Vy (m) is the real and symmetric matrix of expectation values of 
the interaction potential V(r) with respect to our basis functions 4>f\r) for a given £: 

oo 

(m) = [drr 2 V(r) </>f (r) <# (r) . (2) 



In Refs. [[5], [], [7] it has been shown that, for interaction potentials of the power-law form 



V(r) = £ 



a n r 



(with sets of arbitrary real constants a n and b n ), the expectation values (m) can be 
easily worked out algebraically; for their algebraic expression for the most general case, 
consult either Sec. 4 of Ref. jg, or Sec. 3.10 of Ref. f§, or Sec. 2.8.1 of Ref. 0. 

With the aid of the above series expansions, the Fourier-Bessel transformations (|6|), 
and the definition (0) of the matrix elements V^\m) of the interaction potential V(r), 
the matrix Aiij is approximated by the (at least, for all power- law potentials) algebraic 
expression 

N N 

Ma = AK l3 (m) + 2 £ b H (m) V r f{m) + 2m 2 ^ e H {m) V r f{m) 

r=0 r=0 
N N N N N 



+ 2 E E c h VV(m) d S] {m) + m 2 £ £ £ V®{m) e st (m) Vjf\m) 



r(P), 

r=0 s=0 r=0 s=0 t=0 

N N N 



+ E E E C H V sr ] ( m ) Cst V t f\m) . (3) 
r=0 s=0 t=0 



The explicit evaluation of the matrix elements Kij{m) of the kinetic energy squared 
and of the various expansion coefficients 6y(m), cy, and e,ij{m) is a tedious but 

straightforward task. One obtains, for K^m) (cf. Sec. 2.2 of Ref. [0), 

4 m 2 



tt^ + 1)(z + 2) (j + l)(j + 2) 



fc=0 



|r - a| \ r(^ (fc + 1)) T(±(l + r + s+\r-s\- k)) (kj^ 
k J r(|(2 + r + s+ \r-s\)) °° S \ 2 j 



r+^i / . I \ I 7 J. / _L I 1 \ I 7' 1 



r + s + 4 \ T(l(A; + l))r(|(5 + 2r + 2s-A;)) /fc^ 



fc=0 



r(3 + r + s) V 2 
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or, in matrix form, 



K(m) = {Ki,j{m)) = 2m 2 



1 5 
7^ 3 



and, for the expansion coefficients bij(m), 
r / \ 4m 



7T v /(2 + f)( ? + 2) (j + l)(j + 2) 



x EE (-2) 
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=0 s=0 
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\r-s\ , 




_k=0 x 
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fc=0 x 


or, in matrix form, 
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| r _ s | \ £(| (fc + 1)) r(|(2 + r + g + |r - s| - k)) fkji 
k J r(|(3 + r + s+ \r-s\)) C ° S \Y / 
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k ) r(I(7 + 2r + 2 S )) C ° S l 2/ 



b(m) = {bij{m)) 



64 m 

57T 



/ i __l 

3 7\/3 
f If 



7VS 27 
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fn order to evaluate the expansion coefficients cy and ejj(m), we introduce the integrals 



4?V) = 



dk k 2+n 



<pf\k)<pf\k) , n = 0,1,2, 



4rV) 



dfc A; 



2+n 



0l (1) (fc)^ O) (A;) , n = 0,1,2, 



E(k) 



(4) 



the (somewhat lengthy) explicit expressions of these integrals are given in Appendix |B[ 
From the latter, the expansion coefficients cy and e^m) are derived by restricting n to 
the appropriate values: = J^- and e^rn) = lf*\m). Explicitly, these matrices read 



c = (cy) 



f024 
9!!tt 



v/3 



_7_ 
IT 



15 113 y/E 



143 



11 
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e(m) = (e«(m)) 



256 
7!! 7rm 



3^ 



1 

89 
99 



V 



where 



(2n + 1)!! = 1 x 3 x • ■ ■ x (2n- 1) x (2n + 1) , n = 0,1,2,... . 

The analytic results for the expansion coefficients dij(m) have already been derived in 
Ref. |1J; there is no need to duplicate this formula or the matrix d(m) = (dij(m)) here. 

As has been done in Ref. JTJ, we shall adopt, as the simplest model for a confining 
interaction between the bound-state constituents, a linear potential: V(r) = A r, A > 0. 
For this interaction, the general expression forVf (m) given in Refs. 0, ^ simplifies 
to (cf. Sec. 2.2 of Ref. 0) 

\f(m,X) = 



i\ 3 \ 



A 



EE 



\r+s 



r\ si 



r(2£ + i + 3)T(2£ + j + 3) 2m f =0 , „ 

T(2£ + r + s + 4) ; 



i + 2£+2 
i — r 



j + 2£ + 2 
3 -s 



the explicit potential matrices V^(m, A) = (yl-\m, A)) , £ = 0, 1, are given in Ref. |T|]. 

The dependence of all these quantities on the dimensional parameters in the theory, 
viz., the mass m of the bound-state constituents and the slope A of the linear potential, 
may be inferred already on dimensional grounds: 





m) 


= m 2 A%-(l) » 


bij( 


m) 


= mbij(l) , 


dij( 


m) 


= mdij(l) , 




m) 


= ie«(l), 








,A) 


= A V«"(1,1) 
m 



£ = 0,1,... ; 

the expansion coefficients Cij are independent of the bound-state constituents' mass m. 
Factorizing off the dependence on m and A in the matrix M.ij, Eq. (0), we end up with 



N 



Mi 



= Am 2 Kyil) + 2 A ]T [bn(l) + e„(l)] V r f(l, 1) 

r=0 

N N 

+ 2A££c^ 1) (MK i (l) 

r=0 s=0 



\2NNN 

+ — 2 E E E ^(1)^(1.1)^(1)^(1,1; 



r=0 s=0 t=0 
^2 N N N 



+ ^EEE C n ^(1, 1) C st V t f(l, 1) . 
r=0 s=0 i=0 

Approximate solutions of the Bethe-Salpeter equation are found by diagonalizing M. 
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3 A Few Illustrative Results 

Mimicking the analysis of Ref. JIJ, let us investigate first the case d — 1 and N = 0. For 
i — 3 — 0, we find, for the first elements of the matrices K(l), b(l), c, d(l), and e(l), 

„ m „ h m 64 . 1024 .V3 n , 256 

K o(l) = 2, Odo(l) = t=— , c 00 = l — 7^ , Odo(1)=i— , e o(l) 



15 7r 315V3 7T 2 105 n 

and, for the expectation values V"W(1, 1), £ = 0, 1, of the linear potential (cf. Sec. 3 of 



Ref. 



4 0) (l,l) = f, V«(l,l) = ^ 



These matrix elements yield, for the bound-state mass M squared, the analytic result 

„ „ 8896 x 23 / 128 A \ 2 , . . 

M 2 = 8m 2 + -—X + -[- m^O . 

315 7i 7 \4b7Tm J 

From this formula, we obtain, for a bound-state constituents' mass m = 0.9 GeV and a 
typical value A = 0.2 GeV 2 of the slope of the linear potential ||, for the bound-state 
mass M = 2.900 GeV. This is only 10% away from the "exact" result M = 2.637 GeV 
for the ground-state mass, computed for a matrix size d = 15 (that is, a 15 x 15 matrix) 
and N = 49 (that is, taking into account the first 50 basis functions) in the expansions 
performed at intermediate steps. 

In general, that is, for matrix sizes d > 4, the diagonalization of our matrix A4ij can 
be done only numerically. Table [l] illustrates, for the lowest-lying radial excitations, the 
rather rapid convergence of the bound-state masses M, obtained as square roots of the 
eigenvalues of M.ij, with increasing matrix size d to the numerically computed "exact" 
results. 



Table 1: Differences M— 2 m of the eigenvalues M of the instantaneous Bethe-Salpeter 
equation and the sum of the masses m of the bound-state constituents, in units of GeV, 
for two spin-| fermions of mass m = 0.1 GeV, experiencing an interaction described by 
a linear potential with slope A = 0.2 GeV 2 and forming bound states of radial quantum 
number n T — 0, 1, 2 and spin-parity-charge conjugation assignment J PC = 0~ + (called 
l 1 So, 2 1 So, and 3 1 So in the usual spectroscopic notation) as functions of the matrix size 
d, for N = 49 (i.e., considering 50 basis vectors) in the intermediate series expansions. 
The last row compares these differences with the outcome of a numerical computation. 



d 


1% 


2% 


3 1 So 


15 


1.477 


2.147 


2.918 


25 


1.461 


2.095 


2.698 


50 


1.461 


2.074 


2.560 


purely numerical 


1.4613 


2.0740 


2.5580 



Figure |I] shows the dependence of the binding energies (i.e., the differences M—2m 
of the eigenvalues M of the instantaneous Bethe-Salpeter equation and the sum of the 
masses m of the two bound-state constituents) of the lowest-lying bound states on m. 
We observe perfect agreement of our findings with results presented in Fig. 1 of Ref. ||. 



s 



M - 2 m 

3r 

2 . 5 

2 

1 . 5 
1 



0.2 0.4 0.6 0.8 1 "' 

Figure 1: Dependence of the differences M — 2m of the eigenvalues M of the instantaneous 
Bethe-Salpeter equation and the sum of the masses m of the bound-state constituents for the 
three lowest-lying J = ^ bound states, obtained from a time-component Lorentz- vector 
confining interaction kernel involving a linear potential V(r) = \ r with slope A = 0.2 GeV 2 , 
on the bound-state constituents' mass m (all masses in units of GeV). The chosen truncation 
parameters are N = 49 in the intermediate series expansions as well as d = 15 for the matrix 
size (except for m = 0.1 GeV, where, in order to increase the accuracy, d = 50 has been used). 

4 Summary, Conclusions, and Outlook 

In the present investigation we developed a technique for the approximative solution of 
the instantaneous Bethe-Salpeter equation with massive bound-state constituents, by 
a reformulation of this equation of motion as an equivalent matrix eigenvalue problem. 
Combining these findings with the analoguous result obtained within a similar analysis 
for the slightly simpler case of massless bound-state constituents presented in Ref. 
we arrive at the conclusion that, for a suitable choice of basis states in the Hilbert space 
of solutions, it is for a large class of interactions possible to convert the Bethe-Salpeter 
equation in the instantaneous approximation for the involved interaction kernel into an 
eigenvalue problem for an explicitly known matrix, with matrix elements given in form 
of analytic expressions. The main advantage of our formalism is that, due to the scaling 
behaviour of the involved quantities, in actual applications (like fitting procedures) the 
required matrices must be calculated only once (for, e.g., unit values of all the physical 
parameters). As consequence of the explicit knowledge of the matrix representation of 
the instantaneous Bethe-Salpeter equation, the eventual diagonalization of this matrix 
represents the only numerical operation required by this method of solution. 

The next step must be to apply this formalism to realistic (i.e., phenomenologically 
acceptable) models of the interquark forces, capable of reproducing the experimentally 
observed hadron spectra and features by describing hadrons as bound states of quarks. 
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A The "Generalized Laguerre" Basis 

Our choice of basis states for L 2 (R + ) is fixed by the configuration-space representation 



#(r) 



" /exp(-mr)Lf <+2) (2mr) , i = 0,1,2, 



\ T(2£ + i + 3) 



(5) 



which involves the generalized Laguerre polynomials (x) (for the parameter 7) : the 
latter quantities are orthogonal polynomials which are defined by the power series M 



t=o 



1+1 - ,=012 
i-t t\ ' ' ' ' 



and which are orthonormalized, with the weight function x 7 exp(— x), according to 



00 

J dx x 7 exp(— x) (x) (x) 



r( 7 + i + 1) 



i,j = 0,1,2,... 



The necessary normalizability of the Hilbert-space basis functions <pf\r) is guaranteed 
by the positive numerical value of the mass m of the bound-state constituents: m > 0. 
The basis functions (pp (r) defined by Eq. ([5]) satisfy the orthonormalization condition 



drr 2 (j>f ) {r) ( j>f{r) = 8 ij , 1, j = 0, 1, 2, . . . . 



Note that the configuration-space representation of our basis states is chosen to be real. 
The corresponding momentum-space representation (pf\p) of the L2(R + ) basis states 
under consideration is obtained from Eq. (pi) by a Fourier-Bessel transformation (recall 
that one is dealing here exclusively with the radial parts of the /^(-R 3 ) basis functions): 



dppZjtipr)^ (p) , i = 0,1,2,... , £ = 0,1,2, 



(6) 



Explicitly, it reads 

4\p) = 



(-■i)\j=- j drr 2 j e (pr) <j)f\r) , / = 0.1.2. 



£ = 0,1,2, 



(2m) 2e+3 i\ (-i)V 



\ T(2£ + t + 3) 2 i + 1 / 2 v{£ + 
£ (-1)* ( i + 2i + 2 \ Y{2£ + t + 3){2m) t 



t=o 



i-t 



(p2 _|_ m 2\(2£+t+3)/2 
2 \ 



V 



, 2 = 0,1,2, 



x 2 2 ' " ' 2 > 2 + m 2 7 

with the hypergeometric series F, defined, in terms of the gamma function T, by 



(7) 



F(u, V] w; z) 



rw 



E 

n=0 



r(u + n) F(v + n) z r ' 



F(w + n) 



nl 
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The momentum-space basis functions cpf (p) satisfy the orthonormalization condition 

' dpp 2 <pf\p)<j ) f{p) = 5 l3 , z,j=0,l,2,... . 



The availability of the Fourier transform of our basis functions (pf 1 (r) in analytic form 
represents the main advantage of our choice (pj). Note that the momentum-space basis 
functions are real for £ = 0, as well as for all even values of £: 

<pf\p) = ^\p) for £ = 0,2, A,... , Vi = 0,l,2,... . 

In order to get rid of that rather difficult-to-handle hypergeometric series F in Eq. (0), 
we occasionally take advantage of a somewhat simplified form of the momentum-space 
basis functions (j>\ (p), namely, for £ = ||, 

(p) = iE£:E(-2)'(Hi) 



do) 



mnT(i + 3) p t=0 

n \ -(t+2)/2 



i + 2 
i-t 



x 



1 + 



P 



sin ( (t + 2) arctan 



m 



and, for £ 
^(p) 



x 



X 



E 

t=0 



\ tv (i + 1) (i + 2) (z + 3) (i + 4) p 2 
(-2)'/i + 4\ (t + 3)!m* 



j - t I (p 2 + m 2 )(*+ 3 )/ 2 



Vp 2 + rri 2 . 



t + 2 



P 



sin [(t + 2) arctan — 



rn 



m 
t + 3 



sin ( (t + 3) arctan — ] 
V m/ 



The latter form of the basis functions is obtained with the help of a suitable recursion 
formula ||]. 



B Some Useful Integrals 

With the help of the simplified forms of the basis functions <pf^ (p) and (j)^ (p) recalled 
in Appendix [A], explicit expressions for the integrals defined in Eq. may be found: 

4V) 

4m n-l 



7T^/(z + l)(2 + 2) (j + l)(j + 2) 



X 

r=0 s=0 



X 



]r- S ] 

E 

fc=0 



s| \ r(| (k + n + 1)) r(| (4 + r + s + |r - s| - n - k)) 



■ 7T 



COS 



r(|(5 + r + s+ \r- s\)) ^ \ 2 



-+£+ 4 / r + s + 4\ T(l(k + n + 1)) T(| (8 + 2 r + 2 a - n - fc)) 
to I * J " r(|(9 + 2r + 2 S )) C0S {2 / 
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and 



4V) 



i m 
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